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Abstract 

We study variational problems for curves approximated by B-spline curves. We show that one 
can obtain discrete Euler-Lagrange equations for the data describing the approximated curves. 

Our main application is to the occluded curve problem in 2D and 3D. In this case, the aim 
is to find various aesthetically pleasing solutions as opposed to a solution of a physical problem. 
The Lagrangians of interest are invariant under the special Euclidean group action for which 
B-spline approximated curves are well suited. 

Smooth Lagrangians with special Euclidean symmetries involve curvature, torsion, and arc 
length. Expressions in these, in the original coordinates, are highly complex. We show that, 
by contrast, relatively simple discrete Lagrangians offer excellent results for the occluded curve 
problem. The methods we develop for the discrete occluded curve problem are general and can 
be used to solve other discrete variational problems for B-spline curves. 



1. Introduction 

Finding a path which connects two existing paths as described in Figure [1] is a practical 
problem which has applications in physics, engineering, computer vision and many other areas; 
mathematically, the occluded curve problem leads to the study of the Euler-Lagrange equations 
derived from the extremization of a variational problem. 
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As we can see from Figure [TJ there are many possible solutions, even a straight line could 
be a "solution" of the occluded curve. Therefore, certain constrains have to be applied to the 
solution curve: 

- The solution curve has to be smooth, which means the boundary conditions must ensure 
that the occluded curve has C 1 continuity at the boundary; 

- The solution method needs to be equivariant under certain transformation (such as 
rotation and translation), which means that completing the curve before translation and 
rotation is the same as completing the curve after; 

- The solution curve needs to satisfy certain optimization criteria. 

For plane curves, these three conditions can be satisfied simultaneously if the solution curve 
/(s) = (x(s),y(s)) extremizes a functional J" L(s, n, k s ■ ■ ■ )ds, where Lds is the Lagrangian, 
and L is a function of the invariants of the Euclidean group, where s is the arc length, and k is 
the Euclidean curvature of /(s). For space curves, torsion and its derivatives may be included 
in the Lagrangian. 

The solution curve then satisfies the Euler-Lagrange equations for this variational problem. 
The method of solving the occluded curve in the smooth case can be found in the bookpQ. We 
give only a brief review of the method in the next paragraph. 

In the smooth case, the Euler-Lagrange equations are very complex. In 2D, for L — k 2 , the 
solution of the Euler-Lagrange equations which is known as Euler's elastica yields that the 
curvature k is a Jacobean Sine function; in 3D, for the same Lagrangian as in 2D, L = k 2 , 
the solutions of the Euler-Lagrange equations show that the curvature and the torsion are 
Weierstrass p functions. Finding the explicit solutions such as x(s) and y(s) involves yet more 
complicated calculation. So using approximation method becomes important and necessary. 

In computer vision and many other fields, all geometric entities are represented by a variety 
of approximation methods. We use B-spline curves in this paper to approximate given curves 
and develop a discrete variational calculus for the occluded curve problem based on them. 

The input for the occluded curve problem is the outside pieces of the curve /i and fi (see 
Figure [T]). We can use the control points of their B-spline approximations as discrete input data 
for the boundary conditions. Then we can construct the discrete Euler-Lagrange equations, 
which are for the control points of the missing piece of the curve. 

Similarly to the smooth case, the solutions of the discrete Euler-Lagrange equations depend 
highly on the choice of the Lagrangian. Choosing an appropriate Lagrangian, which may or may 
not represents a corresponding smooth Lagrangian, we will be able to obtain an aesthetically 
pleasing result. 

Different approaches of the similar problems are studied by others. Some used discrete points 
on the curve to approximate and construct invariants [8 which has disadvantage when dealing 
with complicated curves, the data need to be processed are massive. Also precision is lost when 
reducing points to simplify the problem. 

The Euler-Lagrange equations for many kinds of discrete variational problems are difference 
equations. Our Euler-Lagrange equations are different. Difference equations are special type 
of recurrence equations which recursively define sequences; our discrete equations, on the 
other hand, define quantities without recurrence. The analysis and calculation of our discrete 
equations are much easier than difference equations and hence grant us advantages in solving 
the discrete occluded curve problem. 

In this paper firstly we give the context of the discrete occluded curve problem. Then a brief 
introduction to B-spline curves and our main tool, syzygies on the control points, which we will 
define. Using discrete syzygies is essential to constructing the Euler-Lagrange equations. After 
that, we introduce the boundary conditions which simplify the discrete system by eliminating 
some unknowns. The simplified system can be then solved. Some examples both in 2D and 3D 
are given, and readers can compare the results from using different Lagrangians and input data. 



DISCRETE VARIATIONAL CALCULUS 



Page 3 of [H 



From those examples we can see that choosing an appropriate Lagrangian is very important of 
obtaining an aesthetically pleasing solution. A general algorithm is then described. 

By considering different sorts of boundary conditions, more features can be introduced into 
the solution curve. Examples are given in 2D space which are classified into different cases 
based on the boundary conditions and input data. Finally, a revised algorithm is given which 
is designed not only for the discrete occluded curve problem but also for the discrete variational 
problem in general. 

Without further ado, let us now introduce the discrete occluded curve problem. 



2. Discrete Occluded Curve Problem 

The discrete occluded curve problem is based on approximation. In the industrial field of 
computer graphic design, using polynomials to approximate actual curves is very practical, and 
B-spline[2] has proved to be one of the most efficient methods due to its outstanding properties. 
We will use B-spline curves to simplify our variational problem especially because of its affine 
invariance property, which implies, that the curve after transformation can be re-constructed 
by the transformed control points. Hence we can convert the differential invariants problem to 
the finite difference invariants problem, which is much easier. 

Talking about invariants, we need to give the group actions which act on our curves. And the 
geometric entities such as curvature and torsion which will be used to analyze the variational 
problem are invariant under those group actions. The group actions we will use in this paper 
are those of the special Euclidean group actions SE(n), which include rotation and translation 
in n dimensional space. 

Notation: I will use • to represent the group action, i.e. x = g • x, g <E SE(n); and ■ to 
represent the inner product between two vectors. 

Considering the input curves are at the "normalized" positions, instead of finding the solution 
curve, we only need to find the control points of the solution curve. This leads us to the problem 
of finding and solving discrete Euler-Lagrange equations. 

2.1. B-spline Curves and Syzygies 

Before we start the discussion of the discrete Euler-Lagrange equations, we will give a brief 
overview of the B-spline basis functions and syzygies pQ. 

In computer graphics and computer aided design area, the B-spline usually refers to the 
parametric spline curves that are constructed by interpolating the control points using the 
B-spline basis functions which are defined over a certain knot sequence. The control points of 
a B-spline curve, just as the name implied, are the points which control the shape of the spline 
curve. The basis functions of a B-spline are defined recursively over a knot sequence, which is 
a sequence of monotone increasing real numbers. When the knots are equidistant, the spline is 
called uniform, otherwise non-uniform. The degree of the spline curve depends on the degree 
of the basis functions. 

Let a B-spline curve be given as 

n 

fj(t) = B(t)P = XX fc (t)pj, t G [t J+k ,t J+k+1 ], j = l..n - k 

where B(t) — (b? k (t)) is a row vector of basis functions which are defined over the knot sequence 
{t\,t2, - ■ ■ ,t m }, P = (pi) is a column vector of control points, j is sequential number of the 
spline curve, and k is the degree of the spline. 
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Notation: When n = k + 1, it is a one-piece spline curve, where j = 1, the subscript of fj (t) 
will be omitted; when n > k + 1, it is a multi-piece spline curve, where j — l..n — k. And m..n 
as in Maple means a set of integers from m to n 

The basis functions are denned recursively over a knot sequence {t\,t2, ■ ■ ■ ,t m }, with < 
U < U + i < 1, Vi, m = n + k + 1, where n is the number of the control points: 



1 ti < t < fj_|_i & tj < fj_|_x 

otherwise 



Let a B-spline curve /(i), i £ [0, 1] with control points (pi)™=i be given. In order to have certain 
properties of the B-spline curve, such as /(0) = pi, /(l) = p n , line segments pipl and p n -\Pn 
are tangent vector at /(0) and /(l) respectively (Bessel end conditions) , we need to set the 
first fc + 1 knots to and the last k + 1 knots to 1 . In the following analysis and examples we 
use the uniform B-spline curves which satisfy the Bessel end conditions. 

Considering the control points of the B-spline curves at the normalized positions as discrete 
invariants, the differences between the control points are difference invariants. Those difference 
invariants will be utilized in the Lagrangian. 

In the variational problem using discrete Lagrangians, in order to find the extremized 
solution, we introduce a dummy variable e, and apply a variation to the difference invariants 
with respect to e. The syzygies are the relations between the variation of different difference 
invariants. 

More definition of syzygies in smooth case can be found in the bookpQ. Following sections 
will give the construction of the syzygies in the discrete case. 



(2.1) 



2.2. Construction of Discrete Euler-Lagrange Equations 

Considering now the control points of the given curves as input data. The discrete Lagrangian 
is defined over the whole interval of control points, but the variation only applies to the control 
points of the occluded curve. The output will be the control points of occluded curves, which 
extremize the variational problem. 

Firstly we discuss the construction of the discrete Euler-Lagrange equations bases on single- 
piece and then multi-piece cubic B-spline curves which satisfy the Bessel end conditions 
respectively. 



2.2.1. Onc-Piece B-spline plane Curve Let a one-piece cubic B-spline plane curve be f(t) = 
Yli=x bi,a(t)Pi, where the basis functions are defined over the knot sequence {0, 0, 0, 0, 1, 1, 1, 1}, 
be a solution to the Euler-Lagrange equations derived from the extremization of the variational 
problem f|2 . 2|) . The corresponding Lagrangian is invariant under the group action g 6 SE(2). 
We have 

C[I,I K ]=Y J L{h,h, K ), (2.2) 

i 

where the invariants I = (7j) = (pi) are the control points at the normalized position, Ik = 
(Ii,k) = (pf) are the differences of control points at the normalized position. Here K = {1,2,3} 
is a set of difference orders. The normalized position is defined as the last control point of left 
input curve locating at the origin, and the first control point of right input curve locating at 
the x-axis (refer to Figure [T]). The difference of the control points are defined as: 

Pi=P l il{-P l f\ Pl=P*, « = l-3 
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Now introduce a dummy variable e and apply a variation to the difference invariants with 
the respect to e, so that If, i = 1.. 4 is the variation of ij, and the syzygies of the difference 
invariants for the control points are: 
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(2.3) 



where S" is the shift map, = ii+i, and id is the identity map. Notice that the shift 

operator and the group action do not commute, which means S(Ii t i) = S(g • (pi+i ~ Pi)) 
g • (S{pi+i — Pi))- The shift operator here only is defined on the normalized difference of the 
control points. Since the variation only happens on the occluded curve, the shift operator and 
its inverse operator will send the variation If, i > n & i < 1 to zero, corresponding to the input 
data being fixed. Hence the extremization of the variational problem (|2.2[) becomes: 
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(2.4) 



dIi,K4 
Boundary Terms) 



If +B.T.S 



The discrete Lagrangian does not need to be exact discretization of smooth Lagrangian. It just 
needs to be invariant and have some geometric meanings, easily solved, and the solution should 
be aesthetically pleasing. Now let us see some examples. 



Example I. Let an occluded curve be one-piece cubic B-spline curve, with the basis 
functions defined over the knot sequence {0, 0, 0, 0, 1, 1, 1, 1} and the control points pi, i = 1..4. 
Choose the Lagrangian L = I\^ • ^2,2, which represents the angle between two second order 
differences of the control points, with I\^ = (S — id) 2 ii and I2.2 = (S — id) 2 /2- Then the 
corresponding Euler-Lagrange equations are: 

E 1 = (S- 1 - id) 2 [(1, 1) T • 7 2)2 ] E 2 = (S- 1 - id) 2 [h a ■ (1, 1) T ] (2.5) 

Example 2. Let an occluded curve be the same as in Example [TJ now let the Lagrangian 
be L = Ji : i • J3 1, which represents the angle between two first order differences of the control 
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points, with Ji i = (S — id) Ji and = (S — id)/3. Then the corresponding Euler-Lagrange 
equations are: 

£ 1 = (S- 1 -id)[(l,l) T -/ 3 ,i] S 2 = (S- 1 -id)[Ji,i-(l,l) T ] (2.6) 

As we can see from the examples above, choosing different Lagrangians will give us different 
Euler-Lagrange equations, thus different solutions to the occluded curve. 

The basis functions of a cubic B-Spline curve are already given, so there are 4 control points 
which contain eight unknowns to determine, but we only got two Euler-Lagrange equations. 
Hence we need some boundary conditions to eliminate the extra unknowns. 

Determine the boundary conditions 



P-l P2 




P-2 

Figure 2. Determine the boundary conditions of the solution path fo 

Figure [5] shows the occluded curve problem with the control points pi = (xi, yi) T ■ As we can 
see the control points (red) of the solution curve must satisfy some boundary conditions: 

1. The first and the last control points (pi and pn) of the occluded curve are the last and 
the first control points of the two input curves (/i and / 2 ) respectively; 

2. The second and the second last control points (p2 and P3) must lie on the tangent line of 
the first and last curve points respectively. 

According to the two rules above, if we set 

P2 = a(pi - Po) P3 = /3(P4 - Ps), (2.7) 

with a, (3 £ R, as shown in Figure^ we can reduce the eight unknowns to two, thus we can use 
the two Euler-Lagrange equations in f|2.5[) and (|2.6[) to solve the control points of the occluded 
curve. 

Since we require the solution curve are calculated at the normalized positon, we can transform 
the original input curves to the normalized position first and calculate the solution curve, then 
transform them all back to the original position. 

Solving the discrete Euler-Lagrange equations 

Solving Example [7] Now let us solve the problem described in Figure Let two 
one-piece cubic B-spline curves f\ and f% been given, their control points are p-2 = 
(0,0),p_i = (l,4),po = (2,1),?! = (4,3),p 4 = (7,2)^5 = (8,3)^ 6 - (9,l),p 7 = (10,2). First 
transform them to the normalized position, that is p\ — (0, 0), 1/4 = 0, after solve the problem, 
transform them back to the original position. Using the Euler-Lagrange equations in \2. 5|) and 
the boundary conditions discussed previously, we will obtain the solution curve as in Figure^ 
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fx 




Figure 3. Solution curve fo to the problem in Figured using Euler-Lagrange equations in 

i)2.5|) is shown in red 

Solving Example [2| Using the same input data and conditions as in Example [TJ but the 
Euler-Lagrange equations in 1^2.6}) . the corresponding solution curve is shown in Figure 



FIGURE 4. Solution curve f to the problem in Figured using Euler-Lagrange equations in 
H'2.6\) is shown in red. The solution curve is smooth, but has very high curvature at the right 

hand end. 

Remark 1. Comparing the two solution curves using different Lagrangians, we can see 
that, the solution in Figure |3] obviously is better than in Figure |U That is because the 
Lagrangian we used in Figure |3] includes the second order difference of the control points 
which relates to the dominant derivatives in the curvature form. However, the Lagrangian we 
used in Figure 0] only contains the first order difference, which is aesthetically less pleasing. 
Hence the choice of the Lagrangian is crucial in solving the occluded curve problem. 

Next We consider the same problem but in three dimensional space. 

2.2.2. Multi-piece Space Curve The construction of the basis functions of a space curve is 
the same as of a plane curve, we only need to add one more dimension to the control points. 

In order to gain more local control and admit more complexity of the occluded curve, 
sometimes we need to use multi-piece B-spline curve to solve the occluded curve problem. 
As for the one-piece plane curve, the key is to find an appropriate Lagrangian. Now let us see 
some examples. 

Example 3. let a two-piece cubic space curve /oi(i) = 2i=i ^1 3^)Pi an d f(aif) = 
2<=2 1% 3 (t)pi be the solution to the occluded curve which connects two one-piece 




h 
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cubic curves /i(t) = J] \ = _ 2 b ii3 (t)pi and / 2 (t) = 2i=5 &*,3(*)Pi, where p_ 2 = (0,0, -2),p_i = 
(1,4, -1.5), po = (2, 1, -l),Pi = (4, 3, -0.5), p 5 = (7, 2, 0),p 6 = (8, 3, 0.5), P7 = (9, 1, l),p 8 = 
(10,2,1.5). The basis functions b} 3 (t) 7 t e [0,1/2] and b 2 l3 (t),t G [1/2,1] are defined over the 
knot sequence {0,0,0,0,1/2,1,1,1,1}. The Lagrangian we chose in this example is L = 
{h,2 x £2,2) • ^3,2 + ^1,3 • ^2,3- The solution is shown in Figure [5] 




Figure 5. Left shows the solution of the occluded curve as a two-piece cubic space curve in 
Example^ red curves show the input curves f\ and fa, blue curve shows foi(t), t e [0, 1/2] 
and green curve shows fo2(t), t 6 [1/2, 1]; right shows the solution of the occluded curve as a 
two-piece cubic space curve in Example^ red curves show the input curves f\ and f%, blue 
curve shows foi(t), t G [0, 1/2] and green curve shows /o2(i)> t 6 [1/2, 1]. 



Example 4. Same as in Example [3j let a two-piece cubic space curve /oi(i) and fo2(t) 
be the solution to the occluded curve. But let the input data now be p-2 = (0,0, 0),p_i = 
(2, 1,1), p Q = (1,2, 1.5), Pl - (0.5, 1.5,2), p 5 = (0,0,3),p 6 = (-2, l,4),p 7 = (1, 2,4.5), p 8 = 
(0.5, 1.5, 5). The Lagrangian we choose in this example is same as in Example |31 The solution 
is shown in Figure [S] 

Now let us generalize an algorithm of solving the occluded curve problem, which is 
approximated by uniform B-spline curves that satisfies the Bessel end conditions: 

Algorithm 1. 

Input: Two sets of control points (jPi)\L ia and ( P j)jLj , the degree of the occluded curve d, 
the number of the segments of the occluded curve I, and a discrete Lagrangian L. 

Output: A set of control points (PkY^ =if which satisfy the extremization of the variational 
problem subjects to the discrete Lagrangian L; 

Step 1: Transform the input curves to the normalized position; 

Step 2: Put the corresponding differences of control points into L, use the method introduced 
in (|2.4p to calculate the Euler-Lagrange equations and obtain the solution; 

Step 3: Transform the input curves and the solution curves from the normalized position 
to the original position. 
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Remark 2. Algorithm [JJ shows us, that for any discrete input data which describe an 
occluded curve problem approximated by B-spline curve, we can find an aesthetically pleasing 
solution curve, which extremizes the variational problem subjects to an appropriate Lagrangian. 

The choice of the Lagrangian is crucial, not only because it yields different results, but it 
also decide the solvability of the discrete system. We need to use the Lagrangian to derive 
sufficient number of Euler-Lagrange equations. For example, using L = l\i ■ I2.2 will derive 
two independent Euler-Lagrange equations, while using L = + -^2,2 will only derive one 
Euler-Lagrange equation. 

The input data we choose are relatively simple, so we only care about the boundary conditions 
described in Figure [5] But what if we have relatively complex input data, will Algorithm [T] still 
offer us satisfied solutions? 



3. Adding Complexities to Solution Curve 

From the construction of the occluded curve we can see, that the complexity of the solution 
curve highly depends on the input data, the degree of the solution curve, and the order of 
difference invariants in Lagrangian. Using one piece cubic spline as solution curve limits its 
complexity. As shown in the previous examples, one piece cubic spline curve will be enough to 
solve the occluded curve problem described in Figure [TJ but what should the solution look like 
in Figure [SP 




FIGURE 6. Adding complexity to a path connecting the paths f\ and $2 



3.f . Complexity of Occluded Curve 

The solution curve to the problem described in Figure |B] requires complexity with more than 
one inflection point. It can be obtained by inserting control points, which will result in either 
increasing the degree of the solution curve and keep the number of spline piece; or increasing 
the number of spline piece and keep the degree of the solution curve. Let us now discuss the 
case of increasing complexity of the solution curve by inserting one control point. 

The complexity of a one piece cubic solution curve will be determined by the input data, 
and can be classified into two cases: 

- Case 1: The tangent vectors at the entering (left) and leaving (right) point of the occluded 
curve have same signs. In this case, the solution curve has one inflection point. An example 
is shown in the left picture of Figure [JJ 

- Case 2: The tangent vectors at the entering (left) and leaving (right) point of the occluded 
curve have different signs. In this case, the solution curve has no inflection point. An 
example is shown in the right picture of Figure [7] 
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Figure 7. Left shows a one piece cubic solution curve with one inflection point, 
corresponding to the same signs of the tangent vectors at the entering and leaving points of 
the occluded curve; right shows a one piece cubic solution curve with no inflection point, 
corresponding to the different signs of the tangent vectors. The Lagrangians used in both 

cases are the same as in Example^ 

3.2. Adding Complexity 

Adding one more control point to a one piece cubic solution curve will either make it become 
a two piece cubic curve or an one piece quartic curve. As we can see in the Figure UJ adding 
one control point to the occluded curve in the left picture will not add inflection point; adding 
one control point to the right picture will add two more inflection points to the occluded curve. 
We need to discuss these two cases respectively. 

3.2.1. Case 1 Let an one piece cubic solution curve be f(t) = J2i=i bi(t)Pi, where the basis 
functions are defined over (0, 0, 0, 0, 1, 1, 1, 1), after inserting one control point, it becomes either 
of the following curves. 

Two piece cubic curve: fj(t) — bj_j +1 (t)qi, j = 1,2, where the basis functions are 

defined over (0,0,0,0, 1/2,1,1,1,1); 

One piece quartic curve: f(t) = Y^=i &i(*)<Z»> where the basis functions are defined over 
(0,0,0,0,0,1,1,1,1,1). 

Using either two piece cubic curve or one piece quartic curve as solution curve we will 
have five control points in total qi — (xi, t/,), i = 1..5, which can construct three second order 
difference invariants. Hence we can use the Lagrangian: L = (Ii t 2 ■ I2.2) {I 2, 2 1 ^3,2), where I^k — 
q^i — q^~ ) are difference invariants at the normalized position, where the last control point 
of the entering (left) curve at the origin and the first control point of leaving (right) curve on 
the a;-axis. 

From the boundary conditions in (|2.7p . we obtain two unknowns a and /3; other two unknowns 
are the coordinates of 93. Since we only have three Euler-Lagrange equations, we need to use 
a boundary condition to eliminate one unknown in q^. Judging from the normalized position, 
noticing that in this case inserting control point q^ will not add inflection point, so we do 
not need to constrain the y-coordinate of (73. Therefore, set X3 = (X2 + X4)/2, solve the three 
Euler-Lagrange equations for a, j3 and j/3, we will obtain the solution curve. Next let us see an 
example. 

Example 5. Input are two cubic curves, one is defined over (0, 0, 0, 0, 1/4, 1/2, 3/4, 1, 1, 1, 1) 
with control points (-6, 0), (-5, -1), (-4, 0), (-3, 1), (-2, 0), (-1, -1), (0, 0); and another is 
defined over (0, 0, 0, 0, 1, 1, 1, 1) with control points (3, 0), (4, 2), (5, 2), (6, 0). A two piece cubic 
solution curve is defined over (0, 0, 0, 0, 1/2, 1, 1, 1, 1) and an one piece quartic solution curve 
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is defined over (0,0,0,0,0,1,1,1,1,1). Using the Lagrangian L = (ii,a • h,2)(h,2 -ha), the 
solution curves are shown in Figure |S] 



Figure 8. The solution curves of Example O Left shows the solution curve as a two piece 
cubic spline curve, whose basis functions are defined over (0, 0, 0, 0, 1/2, 1, 1, 1, 1) at the 
normalized positions; right shows the solution curve as an one piece cubic spline curve, whose 
basis functions are defined over (0, 0, 0, 0, 0, 1, 1, 1, 1, 1) at normalized positions. 



3.2.2. Case 2 All setups are the same as in case 1, except the constrain to eliminate one 
unknown in q^. Since we want to add inflection points to the solution curve, which means 2/3 
need to have different sign as 1/2 and j/4. Set j/3 = — (?/2 + ?/4)/2, solve the three Euler-Lagrange 
equations for a, /3, and X3, we obtain the solution curve. Let us see an example: 

Example 6. Input are two cubic curves, one is defined over (0,0,0,0,1/3,2/3,1,1,1,1) 
with control points (— 2, 6), (— 1, 3), (0, 6), (1, 4), (2, 3), (4, 3); and another is defined over 
(0,0,0,0,1,1,1,1) with control points (7, 2), (8, 1), (9, 3), (10, 2). A two piece cubic solution 
curve is defined over (0, 0, 0, 0, 1/2, 1, 1, 1, 1) and an one piece quartic solution curve is defined 
over (0, 0, 0, 0, 0, 1, 1, 1, 1, 1). Using the same Lagrangian as in Example [5j the solution curves 
are shown in Figure [3] 



Figure 9. The solution curves of Example® Left shows the solution curve as a two piece 
cubic spline curve, whose basis functions are defined over (0, 0, 0, 0, 1/2, 1, 1, 1, 1) at the 
original and normalized positions; right shows the solution curve as an one piece quartic 
spline curve, whose basis functions are defined over (0, 0, 0, 0, 0, 1, 1, 1, 1, 1) at the original and 

normalized positions. 
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3.3. Adding More Constrains 

From the example in case 1 we can see that quartic solution curve looks more smooth than 
the cubic solution curve; from the solution in case 2 we can see, after two inflection points 
are added, the solution curves seem a little too "complex" , the curvature around the inflection 
point is too high. Hence we need to consider the length of the interval between the starting 
and ending points of the solution curve on £-axis. If we increase gap, the solution curve will 
become more smooth, as shown in Figure 1101 

Solution curves of other examples which using the same Lagrangian as in Example [5J using 
one piece quartic and two piece cubic spline curves, with different gaps are shown in Figure [TT] 
and[H 




FIGURE 10. Solution to Example® comparing with Figure® the gap between two input 

curves are increased. 




FIGURE 11. Left shows a two piece cubic solution curve, right shows an one piece quartic 

solution curve 
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Figure 12. Using same input data as in Fisure [7T1 but increase the gap between the two 
given curves. Left shows a two piece cubic solution curve, right shows an one piece quartic 

solution curve 

Assuming that the left input curve has more complexity than the right input curve, now we 
can generalize a revised algorithm: 

Algorithm 2. 

Input: Two sets of control points {pi)\= io and (pj)jLj , the degree of the occluded curve d, 
the number of the segments of the occluded curve Z, and a discrete Lagrangian L. 

Output: A set of control points (Pfe)feLj which satisfy the extremization of the variational 
problem subject to the discrete Lagrangian L; 

Step 1: Transform the input curves to the normalized position; 

Step 2: Use the gap between two input curves as a measure, calculate the number of the 
inflection points n\ and n-z in the same gap length of the left and right input curves which are 
adjacent to the occluded curve. The number of the inflection points for the solution curve is 
n = (ni +rt2)/2, and rounding to the integer part. Depending on the degree of the solution 
curve d, determine the number of spline pieces. Using a one piece cubic curve as basic solution 
curve, adding control points to it. 

Step 3: Judging from different cases introduced previously, using boundary conditions and 
constrains to eliminate some unknowns, put the difference invariants into L. Use the method 
introduced in (12.41) to calculate the Euler-Lagrange equations and obtain the solution; 

Step 4: Transform the input curves and the solution curves from the normalized position 
to the original position. 



We have showed that a complex variational problem in the smooth case can be simplified by 
using an approximation method, such as B-spline curves, and B-spline approximated occluded 
curve can be obtained with relative case and aesthetically pleasing effects. 

When complexity is not required, the solution curve can be approximated by a one piece 
cubic B-spline curve. Higher order B-spline curves are used when more smoothness is required. 

When using a multi-piece B-spline curve as solution curve, we need to find an appropriate 
Lagrangian which can derive a sufficient number of Euler-Lagrange equations to calculate 
the unknowns. For example, a Lagrangian such as L = I\ 2 + -^2,2 will not give us two 
independent Euler-Lagrange equations. When the number of unknowns exceeds the number of 
Euler-Lagrange equations we can obtain, we need to add more boundary conditions. 



4. Conclusion 
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When more complexity is required, a more detailed solution curve can be constructed by 
adding more control points to the original one piece cubic curve. The result will be a multi- 
piece cubic curve or a one piece higher order curve. The additional unknowns can be solved 
for using extra boundary conditions in order to obtain more inflection points in the solution 
curve. 

Transforming the input curves from the original position to a normalized position is not 
necessary, but will help us to simplify the problem. At the normalized position, the complexity 
of the solution curve can be easily analyzed. The normalization process is a simple example 
of a general procedure called "moving frames" , which can be used for general Lie group 
actions. Q] [5] [B] [7]. 

Some extreme cases are not considered in this paper, since the solution are trivial. For 
example, when both input curves have zero slope at the boundary points, the solution is 
obviously a straight line. 

Compared with using a standard numerical solution to a smooth Euclidean invariant 
variational problem [I], the method used here offers more computational advantage, flexibility 
and reliability. 

Further study could focus on the occluded surface problem and discrete variational problems 
other than occluded curve problem. Casting the variational problem in Conformal Geometric 
Algebra will be interesting, challenging, and could lead to further progress. 

This paper is supported by EPSRC grant EP/E001823/1. 
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